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Ph' Abstract 

QJ • We show that grand unified theories based on 50(10) generate naturally the next- 

to-leading baryon number violating operators of dimension seven. These operators, 
which violate {B — L), lead to unconventional decays of the nucleon such as n — )• 
t;;j- , e~ K~^ , e~ TT'^ and p — )• z/vr"*". In two-step breaking schemes of non-supersymmetric 

in 



<S'O(10), nucleon lifetime for decays into these modes is found to be within reach 

of experiments. We also identify supersymmetric scenarios where these decays may 

f^ I be accessible, consistent with gauge coupling unification. Further, we show that 

^^ the {B — L)-asymmetry generated in the decays of GUT scale scalar bosons and/or 

gauge bosons can explain consistently the observed baryon asymmetry of the universe. 

k> ' The induced {B — L) asymmetry is sphaleron-proof, and survives down to the weak 

j^ \ scale without being erased by the electroweak interactions. This mechanism works 

efficiently in a large class of non-SUSY and SUSY S'O(IO) models, with either a 126// 

or a 16i/ Higgs field employed for rank reduction. In minimal models the induced 

baryon asymmetry is tightly connected to the masses of quarks, leptons and neutrinos 

and is found to be compatible with observations. 
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1 Introduction 

Baryon number violation is a very sensitive probe of physics beyond the Standard Model 
(SM). Interactions which violate baryon number [B) are not present in the renormalizable 
part of the SM Lagrangian, but they can arise as effective higher dimensional operators. 
The lowest 5-violating operators [1| have d = 6 and are suppressed by two powers of an 
inverse mass scale. These operators are realized naturally when SM is embedded in a grand 
unified theory (GUT) such as SU{5) and 5*0(10) upon integrating out the heavy vector 
gauge bosons and colored scalar bosons. They lead to the decay of the nucleon into modes 
such as p — 7- e^Tr" and p — ?■ VK^ . Present experimental limits on nucleon lifetime constrain 
the masses of the mediators (vector gauge boson or scalar bosons) to be larger than about 
10^^ GeV, which is close to the unification scale determined from the approximate meeting 
of the three gauge couplings when extrapolated to higher energies. 

An interesting feature of the d = 6 baryon number violating operators is that they 
conserve baryon number minus lepton number {B — L) symmetry, leading to the selection 
rule A{B — L) = for nucleon decay p]. Thus, observation of any decays which violate 
A{B — L) = rule would hint at new dynamics different from those responsible for the 
d = 6 operators. Decay modes in this category include p — )■ z/7r+, n — )■ e~K~^,e~Tr~^ etc, 
which all obey the selection rule A{B — L) = —2. Observation of these decay modes would 
thus furnish evidence against the simple GUT picture with one step breaking to the SM. 

In this paper we study the next-to-leading d = 7 operators, which obey the selection 
rule A{B — L) = —2 for nucleon decay [2], and show that they arise naturally within 5*0(10) 
grand unified theories. In non-supersymmetric 5*0(10) models with an intermediate scale 
we find the nucleon lifetime for decay modes such as n — )■ e~K~^, e~7r+ to be within reach 
of ongoing and proposed experiments. We also identify SUSY 5*0(10) models where these 
decays may be within reach, consistent with gauge coupling unification. While we focus 
mainly on renormalizable 50(10) models with 126^;^ of Higgs bosons employed for rank 
reduction, we show that our results also hold for models with 16^;^ used for this purpose. 

The second main result of this paper is a mechanism for generating the baryon asym- 
metry of the universe at the GUT epoch. The way it comes about is as follows. There 
are heavy scalar bosons and gauge bosons in 5*0(10) theories which generate the d = 7 
operators. These particles have {B — L)-violating two-body decays, which can generate 
the observed baryon asymmetry of the universe naturally, as we show here. This would not 
be possible in the case of {B — L)-preserving decays of GUT scale particles such as the ones 
in 5f/(5). Although grand unified theories were thought to be the natural stage for imple- 
menting the Sakharov's conditions for baryogenesis [3] up until the mid-1980's [1], this idea 
was practically abandoned after the realization that the sphalerons [5] , which violate B + L 
symmetry, would erase any baryon asymmetry that obeyed the A{B — L) = selection rule. 



This is because the effective interactions generated by sphalerons, the non-perturbative con- 
figuration of the weak interactions, are in thermal equilibrium for temperatures in the range 
10^ GeV < T < 10^^ GeV, and violate {B + L) symmetry. However, if baryon asymmetry 
was generated by {B — L)-violating decays of GUT scale particles, they would be immune 
to sphaleron destruction. We show that this mechanism of baryogenesis, which also in- 
duces the d = 7 5-violating operators, is very efficient and occurs quite generically in 
5*0(10) models. In minimal models there is a tight connection between the induced baryon 
asymmetry and the masses of quarks, leptons and the neutrinos. We also note that the 
minimal renormalizable versions of these models [6l[7] have been extremely successful in de- 
scribing neutrino masses and mixings ^, and in particular predicted relatively large value 
for the neutrino mixing angle 6i^, which is consistent with recent results from Daya Bay, 
T2K, Double-Chooz and MINOS experiments ^. The results of the present paper show 
that these models can also explain the observed baryon asymmetry in a manner closely 
connected to the neutrino oscillation parameters. 

This paper is organized as follows. In Sec. 2 we discuss the d = 7 B and (5— L)-violating 
operators. In Sec. 3 we show how these operators arise in unified 50(10) theories, both 
in the non-supersymmetric version and in the SUSY version. In Sec. 4 we address nucleon 
decay lifetime for the {B — L)-violating modes in 5*0(10) models with an intermediate 
scale. Sec. 5 is devoted to GUT scale baryogenesis mechanism tied to the d = 7 operators. 
Here we show the close connection between baryon asymmetry and fermion masses and 
mixings. Finally, we conclude in Sec. 6. 

2 Baryon number violating d = 7 operators 

We begin by recalling the leading baryon number violating operators in the Standard Model 
which have d = 6. There are five such operators with baryon number B = +1 [lll9]: 

0^ = {d'ur{Q,L,)e,„ O2 = (Q.Q,)(«V)*e,„ O^ = {Q^Q,){QkLi)e,,eki 

O4 = {QiQj){QkLi){re),, ■ {Te)M, O5 = (rf'M^)*(nV)* . (1) 

Here we have not shown the color contractions (which is unique in each term via e^/j^), and 
we have suppressed the flavor indices. We have followed the standard notation for fermion 
fields with all fields being left-handed. Thus u'^ stands for the left-handed antiparticle of 
u/j. The spinor indices are contracted via the charge conjugation matrix between fields in 
parentheses. i,j = 1,2 are the SU{2)l indices. The complex conjugate operators of Eq. 
([I]) would of course carry B = —1. 

An interesting feature of the d = 6 baryon number violating operators on Eq. ([1]) is 
that they all carry lepton number L = 1 along with B = 1. Consequently these operators 
preserve (B — L). Thus nucleon decay mediated by these operators would obey the selection 



rule A{B — L) = 0. The decays p -^ e^TT^,VK^ and n -^ V-k^^c^K^ would be allowed 
by this selection rule, while decays such as p — )■ vK^ and n — )■ e~K^, e~7r+, which require 
A{B — L) = —2 would be forbidden. Grand unified theories based on SU{5) and 5*0(10) 
gauge symmetries generate the operators of Eq. ([1]) suppressed by two inverse powers of 
GUT scale masses. With the GUT scale near 10^^ GeV, as suggested by the approximate 
unification of the three gauge couplings, these operators lead to nucleon lifetimes of order 
Xo32 _ xo^^ years for {B — L) conserving modes, which are in the range that is currently 
being probed by experiments. Discovery of nucleon decay into A{B — L) = —2 channels 
such as n — )■ e~K^,e~7T~^ would however suggest that the underlying dynamics is quite 
different from that of the d = Q effective operators of Eq. ([T]). 

As already noted in the introduction, while GUT scale particles can generate a baryon 
asymmetry in their S-violating decays, it was realized that interactions of the electroweak 
sphalerons would wash out any such asymmetry that conserves [B — L). GUT scale baryo- 
genesis thus went out of fashion after the discovery of sphalerons. This was also in part due 
to the leptogenesis mechanism [lOj discovered soon thereafter, which can elegantly explain 
the observed baryon asymmetry with a connection to the small neutrino masses induced 
via the seesaw mechanism |1 Ij . 

Now we turn to the next-to-leading i?-violating operators beyond those of Eq. ([1]), 
which are of dimension seven. These operators are interesting in that they carry {B — L) = 
±2 [2]. While they are suppressed by one additional power of a heavy mass scale, they can 
naturally lead to sphaleron-proof baryogenesis, as we show here. In several instances we 
also find that these operators may lead to observable {B — L) violating nucleon decay. 

There are nine d = 7 baryon number violation operators with B = +1 listed below [21IT2] : 

O, = {d'ur{d'L,yH;e,„ d2 = {d'dr{u'L,yH*e,„ 

O3 = {Q.Q,){d''LkrH;e,,eki, O, = {Q,Q,){d'LkyH;{Te),, ■ ire)^, 

O5 = {Q,e'){d'drH:, de = {d'dr{d'^L,rH„ 

Or = {d'D^dr(La^Q,), ds = {d''D,L,y(d-j^Qi), 

Og = {d'D^dy(d-^^e') . (2) 

We have used the same notation as in Eq. ([1]). Here H is the Standard Model Higgs 
doublet transforming under 5'f/(3)c x SU{2)l x f/(l)y as (1,2, +1/2). D^ stands for the 
covariant derivative with respect to the SU{3)c x SU{2)l x U{1)y gauge symmetry. Note 
that O21 O5 and Oq must be antisymmetric in the down-fiavor indices. Operators of the 
type {d^IpQi){d^Lj)eij are not written, since they are related to those listed in Eq. ([2]) by 
the equations of motion. All vector and tensor operators can be Fierz-transformed into the 
set of operators in Eq. ([2]). 

Note that all operators of Eq. ([2]) carry B = 1 and L = —1, and thus {B — L) = +2, 
with the complex conjugates operators carrying {B — L) = —2. It is these operators which 



can mediate nucleoli decay of the type n — )■ e^K^, e^7r+ and p — > z/7r+. The higher dimen- 
sionahty of these operators would suggest naively that the nucleon decay lifetime would be 
much longer than the ones obtained from Eq. ([1]). However, as we show below, in unified 
theories based on SO{10) with an intermediate scale, these decays may be accessible to 
experiments. Most interestingly, these operators can naturally generate baryon asymme- 
try of the universe at the GUT scale, which is facilitated by the fact that the electroweak 
sphaleron interactions do not wash out a. {B — L) asymmetry generated at such a scale. 

In the supersymmetric version of the standard model, baryon number violation can 
arise through operators in the superpotential analogous to Eq. ([2]). These superpotential 
operators would have dimension six. Holomorphicity of the superpotential would however 
constrain the allowed operators. There is a single operator of dimension six given by the 
superpotential coupling 

O^susY ^ d'd%'L,{Hu),e,„ (3) 

which carries B = —1 and B — L = —2. This operator must be antisymmetric in the down- 
flavor indices owing to Bose symmetry. H^ here is the up-type Higgs doublet of MSSM. In 
the superpotential, this operator will appear with two inverse powers of a heavy mass scale. 
Terms in the Lagrangian resulting from Eq. ([3]) would have two fermion fields, one Higgs 
field and two superpartner scalar fields, for example. When the superpartner scalar fields 
are converted to standard model fermions by a gaugino loop, effective d = 7 operators of 
Eq. (E]) would be generated, suppressed by a factor (M^Msusy)"^, rather than M~^ that 
occurs for Eq. (|2]) without SUSY, where M is the heavy mass scale. Therefore, potentially 
these SUSY contributions can be more significant for nucleon decay. 

We can now present the complete list of {B — L) = —2 effective operators through d = 7 
in the SM by adding to Eq. ([2]) operators with B = 0,L = 2. These operators have been 
classified in Ref. [13]. While not directly related to nucleon decay, these operators arise 
along with the d = 7 operators of Eq. ([2]) in 5*0(10) unified theories, and they are also 
relevant for GUT scale baryogenesis. Here we collect the linearly independent set of these 
operators through d = 7. The leading operator of course is the well-known d = 5 seesaw 
operator [1] 

Od=5 = {ULj) HkHi (re),, ■ (fe)^ . (4) 

The next-to-leading operators are oi d = 7, and there are ten of them, as listed below. 

0[ = {L,L,){Lke') Hi e^jeM, O', = {ULj){L,,e^) Hi {ft),, ■ (re),;, 

O', = {L,L,){Qkd') Hi e^jtM, O', = {ULj){Q,,d') Hi (re),, ■ (re),;, 

O'^ = {LiQj){L,d') Hi eijeM, O', = {UQ,){Lud') Hi (re),, ■ (re),,, 

O', = {L,L,){Qkuy Hk e,„ O'^ = {L,Lj){Qkuy Hi (fe)„ ■ (re),,. 



We shall see the appearance of some of these operators in the embedding of Eq. ([2]) in 
5*0(10) models. In the supersymmetric standard model, the four holomorphic operators 
0[ — O'^ would be allowed in the superpotential (with no significance attributed to the 
spinor contractions of Eq. ([5]) when applied to the superfields) . 

3 Origin of d = 7 S— violating operators in SO (10) 

In this section, we show that the d = 7 baryon number violating operators of Eq. (jj]) 
arise naturally in the context of 5*0(10) unified theories after the spontaneous breaking 
of {B — L), which is a part of the gauge symmetry. The {B — L) symmetry may break 
at the GUT scale so that 50(10) breaks directly to the Standard Model gauge symmetry, 
or it may break at an intermediate scale Mj below the GUT scale. In the latter case 
the intermediate symmetry could be one among several possibilities: SU{A)c x SU(2)l x 
SU{2)r- SU{A)c X SU{2)l x U{1)r; SU{3)c x SU{2)l x SU{2)r x U{1)b^l; SU{3)c x 
SU{2)l X U{1)r X U{1)b-l', or 5*f/(5) x t/(l), with or without left-right parity symmetry. 
In the non-supersymmetric version an intermediate scale is necessary to be compatible with 
gauge coupling unification |il4j, while with supersymmetry the direct breaking of 50(10) 
down to the MSSM is preferable. Even in the latter case, there is room for intermediate 
scale particles, provided that they form complete multiplets of the 5?7(5) subgroup, since 
such particles do not spoil the unification of gauge couplings observed with the MSSM 
spectrum. 

To see how the d = 7 operators of Eq. ([2]) arise within 5*0(10), we first focus on the 
scalar-mediated operators and write down the Yukawa couplings in the most general setup. 
The Higgs fields which can couple to the fermion bi-linears 16jl6j are IOh, 126// and 120//, 
with the couplings of the 10// and 126// being symmetric in flavor indices {i,j) and those 
of the 120// being antisymmetric. The terms in these Yukawa couplings that are relevant 
to the generation of the d = 7 operators are given below [T5|[T6]. 



£(16^6,10//) = h. 



V 



IQj + u'rL,) h - {dlQ, + elLj) h + {^y^iQj + ule) - d1 



+ {mid] + Q,L^) o;^] , (6) 

£(16,16,126//) = h, [{ulQ, - 3u^Lj) h - {d^Q, - Se^Lj) h 

+V3^ {^QiQj - Uie^j + T^-d"^ wi + V3i{QiLj - tu^d'^ ool 

+ VQid'ru] + u're'j) U2 + 2V3i d^ Lj p - 2^fli vl Qj p + 2^3 < v] r] 

-2VSi < Lj X + 2V3t el Qj x - 2^3 d^ e" 6 + VQiQi Lj $ + ....] ,(7) 



£(16^16,120^) = g^J {diQ^ + e'rLj)h-{u'rQj + u^Lj)h,-V2QiLjco1 



I , ,.^ „ .. X ^ . I 



-i e rf.'c?" ri + 2i u'ru^ ri + V2ie d^w" w^ + ^2 i {d'ru^ - e\u^^ Wa 



— ^-^QiQ^^ - ^QiLj^ - 2i dle]b ^ i eu\u]b ^ 
v2 



(8) 



These terms are written in terms of the Standard Model decomposition of the sub-multiplets. 
We have followed the phase convention of Ref. [I5]. e stands for the S'f/(3)c tensor ea^^. 
Here we have not displayed terms that are irrelevant for inducing the d = 1 baryon num- 
ber violating operators. (Specifically, we have omitted color singlet, color octet, and color 
sextet couplings.) The Yukawa couplings obey /ijj = /i^j, /jj = /jj and gi^ = —Qji- The 
S'f/(3)c X SU{2)l X U{1)y quantum numbers of the various sub-multiplets are given as 
follows. 

/i(l, 2, +1/2), Ml, 2, -1/2), a;(3, 1, -1/3), uj\?>, 1, 1/3), 

p(3,2,l/6), p(3,2,-l/6), r/(3, 1,2/3), rj(3, 1, -2/3), 

$(3,3,-1/3), $(3,3,1/3), x(3,2,7/6), x(3,2,-7/6), 

5(3,1,-4/3), 5(3,1,4/3). (9) 

Different fields with the same SM quantum numbers appear in some couplings, they are 
distinguished by subscripts 1, 2 etc. We have used the same notation for fields with the same 



SM quantum numbers in 10^, 126// and 120//, but it should be understood that these are 
distinct fields. After GUT symmetry breaking various subfields with the same SM quantum 
number would mix. Some of these mixings would involve the vacuum expectation value 
of the SM singlet field from the 126//, denoted by A^ carrying (B — L) = —2. It is this 
field that supplies large Majorana mass for the right-handed neutrino through the coupling 
fij^^vlv'^-IS.'^. With (A*^) 7^ 0, trilinear scalar couplings of the type p*ujH, r]*pH, p*^H and 
X*TiH, will develop. This issue will be addressed in more detail below, but we note that 
such couplings are invariant under the unbroken SM gauge symmetry. When combined 
with the Yukawa couphngs of Eqs. (I6])-(|8]), they would induce the d = 7 baryon number 
violating operators of Eq. ([2]) . 

To see how the d = 7 operators arise in more detail, let us focus of the flavor symmetric 
Yukawa couplings of Eqs. (|6])-(l7]). These couplings generate two of the d = 7 operators as 
shown in Fig. [1] Here H{1, 2, 1/2) is the SM Higgs doublet, which is a linear combination 
of the h{l, 2, 1/2) and h (1, 2, —1/2) fields from 10//, 126// as well as any other Higgs sub- 
multiplet with the quantum number of (1, 2, 1/2) in the theory with which these fields 
mix. Similarly, a;(3, 1, —1/3) generically stands for any linear combination of u and {oo'^)* 
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Figure 1: Effective baryon number violating d = 7 operators induced by the symmetric 
Yukawa couplings of 10^ and 126^ of 5*0(10). Here the SM quantum numbers of the 
various fields are a;(3, 1, —1/3), p(3, 2, 1/6), and H{1, 2, 1/2). 



from the IOh, 1^1,1^2 and (wi)* from the 126//, etc. Before estimating the strength of these 
operators, let us examine the origin of the trilinear scalar couplings that appear in these 
diagrams in 5*0(10). 

To see the origin of p*uH and similar vertices, let us recall first the decomposition of 
various SO(IO) fields under the subgroups SU{A)c x SU{2)l x SU{2)r and SU{5) x f/(l). 
Under G(2,2,4) = SU{2)l x SU{2)n x 5f/(4)c, we have the following decomposition: 

16 = (2, 1,4) + (1,2, 4) 

10 = (2, 2,1) + (1,1, 6) 

126 = (1,1, 6) + (3, 1,10) + (1,3, To) + (2, 2, 15) 

120 = (2,2,1) + (1,1,10) + (1,1,T0) + (3,1,6) + (1,3,6) + (2,2,15) . 



(10) 



Under the G(5, 1) = SU{5) x f/(l) subgroup various fields decompose as follows. 

16 = l(-5) + 5(3) + 10(-l) 

10 = 5(2) + 5(-2) 

126 = l(-10) + 5(-2) + 10(-6)+T5(6)+45(2) + 50(-2) 

120 = 5(2) + 5(-2) + 10(-6)+l0(6)+ 45(2) + 45(-2) 



111 



Now, the quartic coupling (126)^, which is invariant (there is a single such coupling), 
contains the term (2, 2, 15) ■ (2, 2, 15) • (1, 1, 6) ■ (1, 3, TO) under G(2, 2, 4). The p*(3, 2, -1/6) 
field is a subset of (2,2,15) fragment. The H{l,2,l/2) field is part of (2,2,15), while 
ci;(3, 1, —1/3) C (1, 1, 6) and one such field is also part of (1, 3, 10). Thus one sees that the 
coupling (2, 2, 15) ■ (2, 2, 15) ■ (1, 1, 6) ■ (1, 3,l0) would contain the term p*ujHA^, where A^ 
denotes the SM singlet field from 126h that acquires a GUT scale VEV. In the G(5, 1) 
decomposition, p*(3,2,-l/6) C T5(6), if(l,2,l/2) C (45,2), and a;(3, 1, -1/3) C (45,2). 
Thus (126)^ contains the term T5(6) ■ (45, 2) ■ (45, 2) ■ I(-IO), which has the piece p*ujH~^. 
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There are three non-trivial invariants of he type (126)^ ■ (126*)^. These couphngs also 
contain the term (2, 2, 15) ■ (2, 2, 15) ■ (1, 3, 10) ■ (1, 1, 6), as can be seen by examining the 
decomposition under G(2,2,4) and separately under (j(5, 1)|j In an analogous fashion one 
sees that the coupling (126)^ ■ (126*10) contains (2, 2, 15) • (2, 2, 15) ■ (1, 3, 10) ■ (1, 1, 6), with 
the (1, 1, 6) arising from the 10. To complete this discussion we also note that there are three 
invariants of the type (120)^ ■ (126)^ which contain (1, 3, 6) ■ (1, 3, 10) ■ (2, 2, 15) ■ (2, 2, 15), 
with the (1, 3, 6) and one (2, 2, 15) taken from the 120. The (1, 3, 6) fragment contains 
the ?7(3, 1,2/3) of Eq. (Q, which would enter the d = 7 operators arising by integrating 
out the flavor antisymmetric 120 fragments. In particular, this term would induce the 
needed p*rjH* vertex. The (1, 3, 6) ■ (1, 3, 10) ■ (2, 2, 15) ■ (2, 2, 15) term also contains x*vH 
vertex, with x C (2,2,15), r] C (1,3,6) from the 120. The (120)^ ■ (126)^ invariant also 
contains the term (3, 1, 6) ■ (1, 3, 10) ■ (2, 2, 15) ■ (2, 2, 15), with the (3, 1, 6) and one (2, 2, 15) 
taken from the 120. This piece generates the terms ^*pH* vertex, with $* C (3,1,6). 
Finally, there are two invariants of the type (120)^ ■ 126 ■ 126*, which also contain the term 
(1,3, 6) -(1,3, 10) -(2, 2, 15) -(2, 2, 15). 

In order to complete 5*0(10) symmetry breaking, additional Higgs fields such as a 45, 
54 or a 210 is needed. The interactions of these fields with 126 and 120 can generate 
further trilinear and quadrilinear couplings. Take for example the case of a 54 employed 
for completing the symmetry breaking. The trilinear couplings (126)^54 and (126)^54 are 
then invariant. Noting that under SU{5) x f/(l) subgroup, 54 = 15(4) + T5(-4) + 24(0), 
we see that the latter coupling would contain a term 1(10) ■ 15(— 6) • 15(— 4). This has 
a piece A'^(l, 1, 0)y26 p(3, 2, 1/6)y26 P*(35 2, —1/6)54, which would mix the two p fields once 
(A*^) 7^ develops, breaking (B—L) by two units. (Note that the p^ge ^^^ P54 carry different 
{B — L) charges.) The cubic coupling (126)^ 54 also contains the terms 5(2) ■ 5(2) ■ 15(— 4) 
and 45(— 2) ■ 45(— 2) ■ 15(4) under SU{5) x f/(l). These terms contain the couplings 
a;(3,l,-l/3)T26 /^(l, 2,1/2)126 P*(3, 2, -1/6)54 and u^ h{l,2,-l/2)^ p(3, 2, 1/6)54 re- 
spectively. These are the desired trilinear couplings for the generation of the d = 7 opera- 
tors. In this case, the mixing of p*(3, 2, —1/6)54 with the p(3, 2, 1/6)^26 i^ utilized in order 
to connect the p field with the fermion fields in Fig. [1] Similar results follow from the 
quartic couplings (126)^ (54)^ and (126 ■ T26) (54)^. 

While trilinear scalar couplings of the type r]*pH and ^pH* do arise for the 126// sub- 
multiplets, these couplings do not directly lead to baryon number violation. The 77 field 
from 126// has the coupling rju'^v^, while $ has the coupling ^QL (see Eq. ([7])). The 



exchange oi rj — p from 126// would lead to an effective operator {u^v^){d'^Li)H.jeij^ while 
that of $ — p would generate the operator {QiLj){d'^Lk)Hi{fe)ij ■ {fe)kh which is operator 



iThe (126)2(126)2 coupling would also contain a term (2,2,15) • (2,2,15) • (1,3,10) • (1,3,10) under 
G(2,2,4), which has a p*ujHA'^ term in it. However, the w appearing here is from the (1,3, 10), which is 
LJ2 of Eq. ([7]) . UJ2 coupling by itself does not violate baryon number |17| , as is evident from Eq. ([7]) . 
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Figure 2: Effective d = 7 baryon number violating operators obtained by integrat- 
ing our fields from the 120. Here the SM quantum numbers of the scalar fields are: 
H{1, 2, +1/2), u{3, 1, -1/3), p(3, 2, 1/6), r]{3, 1, 2/3), $(3, 3, -1/3), and x(3, 2, 7/6). 



Cg of Eq. ([5]). The r7*p_f/' and ^pH* couplings of the 126h sub-multiplets would however 
be relevant for GUT scale baryogenesis, since they also violate {B — L) symmetry. 

In Fig. [2] we display the effective d = 7 operators obtained by integrating out the 
flavor antisymmetric 120 coupling to fermions. We see that four operators are induced this 
way. As already noted, the required trilinear vertices to complete these diagrams arise from 
quartic couplings. It is also possible to replace the p fields in Fig. [2]by a p field from the 



126, in which case the sum of such diagrams would have no definite symmetry property in 
two of the flavor indices. 

Note that all of the d = 7 operators arising from Fig. [1] and [2] respect B — L symmetry, 
as can be seen by assigning {B — L)(A'^) = —2. 

It should be noted that operators O-j — Og of Eq. (|2]) do not arise at tree level by 
integrating out superheavy particles. They can arise via loops, with suppressed strength. In 
the estimate of nucleon lifetime these operators play a subleading role, since the amplitude 
for the decay would be further suppressed by the nucleon momentum (rather than the 
electroweak VEV for operators Oi — Oq). 

What if the Higgs field employed for reducing the rank of 5*0(10) is a 16h rather than 
a 126j:/? In this case the d = 7 diagrams of Fig. [1] and [2] would still arise, albeit in a slightly 
different way. The {B — L) quantum numbers of the Higgs doublets h{l, 2, 1/2) from the 
IOh, 126h and 120h fields are all zero. The IQh contains a SM singlet filed with B — L = +1 
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which acquires a GUT scale VEV. It also contains a h{l, 2, —1/2) field with B — L = —1. 
Similarly, 16h contains a SM singlet field with {B — L) = —1 and a /i(l,2, 1/2) field with 
{B — L) = +1. The trilinear scalar couplings of the type IGh^Qh^^h and 16^16^10^ 
would mix the B — L = Higgs doublet h(l, 2, 1/2) from the 10^ and the h{l, 2, 1/2) Higgs 
from the 16h which has B — L = +1. The light SM Higgs doublet then would have no 
definite B — L quantum number. The (1,2,4) component of 16h under G(2,2,4) contains 
the field p*(3,2, -1/6), and the (2, 1,4) of l6^ under G(2,2,4) contains uj{3, 1, -1/3), and 
thus the coupling p*uH is generated via the IGh^Gh^^h coupling. One could also take 
w(3, 1, -1/3) from the 10^, p*(3, 2, -1/6) from the TGh and H{1, 2, 1/2) from the 16^:^ to 
generate the p*ujH coupling from 16/^-16^10^. It should be noted that the p*(3,2, —1/6) 
field from the 16^ does have Yukawa couplings to fermions, since in this type of models the 
heavy Majorana masses of the v^ fields arise from the couplings 16il6jl6Hl6//, and upon 
insertion of one VEV for the SM singlet field here, the coupling of p* to fermions is realized. 
We point out that the p(3, 2, 1/6) field is partly in the Goldstone mode, associated 
with the breaking of S'O(IO) to SU{^). However, since other fields such as 45, 54 or 210 
should be employed to complete the symmetry breaking down to the SM, one such physical 
p(3, 2, 1/6) will remain in the spectrum. This is because the 45, 54 and 210 all contain a 
p(3, 2, 1/6) field, and only one p is absorbed by the gauge multiplet. 





(c) (d) 

Figure 3: Origin of c/ = 7 baryon number violating operators via the exchange of vector 
gauge bosons 1/q(3, 2, 1/6) and K-(3, 1, -2/3) of 50(10). 

The d = 7 operators of Eq. (^ can also arise by integrating out the gauge bosons 
of 5*0(10). These vector operators are related by Fierz identities to the scalar operators 
displayed in Eq. (|2]). The relevant diagrams are shown in Fig. [31 (Fig. |3] (d) conserves 
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B, but violates L and B — L [13]. The effective operator from this diagram, after a Fierz 
rearrangement can be identified with O'^ of Eq. ([5]).) In these diagrams Vq denotes the 
gauge boson with the SM quantum numbers (3,2,1/6) (same as the quark doublet Q), 
and Vu<= is the gauge boson that transforms as (3, 1, —2/3), the same as that of u'^ fermion. 
In Fig. [3], each vertex conserves i? — L as it should, which can be seen by assigning 
{B — L) = (—2/3, —4/3) to (Vq, V^c). The covariant derivative for the rank reducing field 
126h would contain the term VgVu'^HlA'^y which enters the d = 7 operators. When 16^ 
is used instead of the 126^, the covariant derivative would contain a similar term, but now 
the H and A'^ fields would carry {B — L) = +1 and —1 respectively. 

3.1 S— violating d = 6 superpotential operators in SUSY SO (10) 



d'' 
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T] 



d" 
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u- 



Figure 4: Superfield diagrams for the generation oi d = Q baryon number violating operator 
of Eq. ([3]) in S'O(IO). Diagram (a) arises from integrating out colored fields from 10^ and 
126/f , while diagram (b) arises from 120//. 



We have identified in Eq. ([3]) a single holomorphic operator of dimension six in the 
superpotential which violates B and {B — L). In Fig. |l]we show how this operator can 
arise in SUSY 5*0(10). Fig. H] (a) is obtained by integrating out colored fields in the IOh 



and 126 H, while Fig. H] (b) is obtained by integrating out such fields from the 120 h- These 
are superfield diagrams, effectively generating F-terms in the Lagrangian. Note that the 
effective superpotential is antisymmetric in down-quark flavor in both diagrams (a) and 
(b). This can be explicitly verified, by making use of the eu'^d'^u'^ vertex in Fig. H] (a) 
which is color antisymmetric. 

In Fig. m (a), the vertex Ld'^p is contained in the Yukawa coupling of Eq. ([7]) and the 
vertex u'^d'^u'^ can be either from Eq. ([6]) or from Eq. ([7]). The u'^u and fip transitions occur 
through [IOhY ("^^ 126// 126//) mass term, for example. As for the vertex upH^ here is 
a concrete example. The coupling 126// 10// 210// is invariant, and contains such a term, 
where u C 10//, p C 126// and H^ C 210//. Note that H^ C 210// carries {B - L) = +2, 
since it is part of the fragment (2,2,10) under (7(2,2,4). Of course, the MSSM field H^ 
in this case is an admixture of h{l,2, 1/2) from 10//, 126//, 126// and 210// and carries no 
definite {B — L) charge. 
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As for Fig. 111(b), the vertices xLu'^ and rjd'^d'^ are contained in the 120^ Yukawa cou- 
phng of Eq. (^. The vertex rfxHu can arise from the superpotential couphng (120/^)^ 210/^, 
which is invariant. 

Suporpotential operators of Eq. (|3]) can also arise in SUSY S'O(IO) models which utilize 
low dimensional Higgs representations, for e.g., {10// + 10'^^ + 16^ + 16^^ + 16^ + 16^:^+45//}. 
Such models have been widely discussed in the literature [T81421] . While -R-parity is no 
longer automatic in these models (unlike in the case of 126j/ models), it can be ensured 
by a discrete Z2 symmetry which distinguishes 16^ from the chiral fermions 16, and which 
remains unbroken. The doublet-triplet splitting problem can be addressed without fine- 
tuning in these models via the Dimopoulos-Wilczek mechanism [22]. Generation of fermion 
masses in these models relies on higher dimensional operators such as 16il6jl6Hl6// which 
induce heavy Majorana masses for the v^, and 16jl6jl6//16^ which induce lighter family 
masses and CKM mixings. (16^ and IQh acquire GUT scale VEVs along their SM singlet 
components, while 16'^:^ and 16^^^ do not. The SU{2)l doublet components of 16h and 16^ 
acquire weak scale VEVs, see for e.g., discussions in Ref. [21].) To see the origin of Eq. 
(jS]) in such models, consider generation of Fig. S] (a). The vertex Ld'^ p arises from the 
Yukawa coupling 16j 16^16^1 6^ after a GUT scale VEV is inserted for the 16 h (note that 
16^ contains a p field), and the vertex u'^ d'^ u^ is contained in the coupling 16jl6jl0j/. 
Now, u'^ can convert itself into u C 10'^:^ via the coupling 10// 10'^:^ 45// when the {B — L)- 
preserving VEV of 45// is inserted, while the p can transition into p C 45// via the coupling 
1&^ 16// 45// with the insertion of a 16// VEV. The pa; H^ term is contained in the coupling 
10// 10^ 45//, which completes the diagram. Explicit models [21] contain all these terms 
necessary for generating the rf = 6 superpotential operator. 

4 {B — L)— violating nucleon decay rates in S'O(IO) 

Before discussing the rates for {B — L)-violating nucleon decay, let us note that cer- 
tain scalar bosons and certain gauge bosons would induce the more dominant d = Q 
baryon number violating operators of Eq. ([T]). Scalar bosons with SM quantum numbers 
c<;(3, 1, —1/3), $(3,3, —1/3) and (5(3, 1, —4/3) and vector gauge bosons with SM quantum 
numbers X(3,2,— 5/6) and Vq(3,2, 1/6) can can induce these {B — L) preserving rf = 6 
operators. Current nucleon lifetime limits restrict the masses of these gauge bosons to be 
larger than about 10^^ GeV, and those for the scalar bosons to be heavier than about 10^^ 
GeV (for Yukawa couplings of order 10^^). Since some of the d = 7 operators arise in 
5*0(10) via the exchange of these particles, the above limits have to be met in our estimate 
oi d = 7 decay rates. 

We proceed to estimate the lifetime of the nucleon for its {B — L)-violating decays in 
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so (10) a straightforward wayo The diagrams of Fig. [T]lead to the estimate 



r(n —7- e" 
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Here we have defined the Yukawa couphngs of uj and p fields appearing in Fig. [T]to be Yqq^ 
and Yid-^p. These couphngs are hnear combinations of the Yukawa couphng matrices /i, / 
and g of Eqs. (jH])-® with flavor indices corresponding to the first family fermions. The 
factors D and F are chiral Lagrangian factors, D c:^ 0.8 and F ~ 0.47. I3h — 0.012 GeV^ 
is the nucleon decay matrix element [21], vr = (A^), and v = {H^) = 174 GeV. We have 
defined the trilinear coupling of Fig. [1] to have a coefficient Xvr. As expected, the rate is 
suppressed by six powers of inverse mass, owing to the higher dimensionality of the effective 
operator. The mass of ci;(3, 1, —1/3) is constrained to be relatively large, as it mediates d = 6 
nucleon decay. For Y ^ 10~^, M^^ > 10^^ GeV must be met from the d = 6 decays. As an 



illustration, choose Yqq^ = Y^d'^p = 10 , M^ = 10 
GeV in Eq. (IT^ . This choice would result in r, 
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GeV, M„ = 10^ GeV, and Xvr = 10 
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3 X 10 yrs. Such a spectrum is 



motivated by the intermediate symmetry G(2, 2, 4) = SU{2)l x SU{2)r x SU{4)c (without 
discrete Parity), which is found to be realized at Mj ^ 10^^ GeV from gauge coupling 
unification [H]. As a second example, take Mp = 10^ GeV, M^ = 10^^ GeV, Xvr = 10^^ 



GeV, Yqq^ = Yid'^p = 3 X 10 . This choice of spectrum leads to r„ ~ 4 x 10 yrs 




6 9 12 

Logio[n/GeV] 




6 9 12 

Logio[n/GeV] 



Figure 5: Unification of the three SM gauge couplings obtained with a light p(3, 2, 1/6) and 
two (1,3,0) scalar multiplets at M/ = 10^ GeV (right panel). The left panel corresponds 
to having one p(3, 2, 1/6) and one (1, 3, 0) scalar fields at the weak scale. 

This second choice for the spectrum can be motivated as follows. The unification of 
gauge couplings may occur without any particular intermediate symmetry, but with certain 
particles surviving to an intermediate scale. We have found two examples of this type 
where the d = 7 nucleon decay is within observable range. Suppose the p(3, 2, 1/6) particle. 



"^{B -\- L)-preserving nucleon decay has been studied in the context of i?-parity breaking SUSY in 
Ref. EH. 
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along with a pair of (1,3,0) scalar particles (contained in the 45^, ^Ah or 210^ needed 
for symmetry breaking) survive down to Mj = 10^ GeV. The three SM gauge couplings 
are found to unify at a scale Mx ~ 10^^ GeV in this case. This is shown in Fig. |5] 
right panel. This figure is obtained by the one-loop renormalization group evolution of 
the SM gauge couplings above Mj with beta function coefficients bi = (62/15, —2, —20/3), 
where dgi/dlnfi = bigf / {16tt'^) . Since the d = 6 nucleon decay lifetime would also be 
near r(p — )■ e^vr^) ^ 10^^ yrs with such a unification scale, this scenario would predict 
observable rates for both the {B — L)-conserving and {B — L)-violating nucleon decay 
modes. Another possible scenario for consistent gauge coupling unification is to assume 
that p(3,2, 1/6) and one (1,3,0) scalar multiplet survive down to the weak scale. Again, 
the gauge couplings unify around Mx ~ 10^^ GeV, as shown in the left panel of Fig. [51 
Here the /3-function coefficients used are hi = (62/15,-7/3,-20/3), appropriate for this 
spectrum. The estimate r(?T, — > e~K'^) ^ 10^^ yrs would follow, as in the second example 
above, but now with Yqq^ = Yi^'^p = 10""^. 

Nucleon decay rates arising from diagrams of Fig. |2] are similar, but with a signifi- 
cant difference. Take for example Fig. [2] (b) with the rip*H trilinear vertex. Neither the 
7/(3, 1, 2/3) nor the p(3, 2, 1/6) field would mediate d = 6 nucleon decay, and may be con- 
siderably lighter than the GUT scale. In non-supersymmetric S'O(IO) models it is natural 
that some scalars survive down to an intermediate scale. Typically the intermediate scale 
is of order 10^° — 10^^ GeV if the associated symmetry contains SU{4:)c- Thus it is possible 
that both 1] and p have masses of order 10^° GeV. Nucleon decay rate is now estimated to 

Here Af/j is defined as the coefficient of the trilinear scalar vertex. For Mp = M^ = 10^*^ 
GeV, vn = 10^^ GeV, Y = 10^^, r„ f^ 3 x 10^^ yrs, which is in the observable range. 

The gauge boson exchange diagrams of Fig. [3] would also induce {B — L)-violating 
nucleon decay. However, we find the rates for these decays to be suppressed. Take Fig. 
|3] (a) for example. The vector gauge boson Vq must have mass of order 10^^ GeV, since 
it mediates d = 6 nucleon decay, while V^^c could be much lighter as it does not have 
5-violating interaction by itself. The VqVucH vertex necessary for connecting the d = 7 
diagram has a coefficient of order q'^vr, while the mass of V^c is of order gvR. The amplitude 
for d = 7 nucleon decay arising from Fig. [3] is then given hj A ^ [q'^vrv) / {My My ^) k, 
{g'^/My ){v/vr). This amplitude is a factor {v/vr) smaller compared to the standard d = Q 
nucleon decay amplitude originating from GUT scale gauge bosons. Unless the {B — L)- 
breaking scale vr is close to the weak scale f , not a likely scenario based on gauge coupling 
unification, nucleon lifetime for {B — L)-violating modes from these diagrams would be 
beyond the reach of experiments. However, as we shall see in the next section, the {B — L)- 
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violating decays of the gauge boson can naturally explain the observed baryon asymmetry 
of the universe. 

As for the supersymmetric diagram of Fig. |H the d = 6 operator in the superpotential 
Eq. (|3]) has a strength A = {Yucd'^^i'^YLdcpY^-pH^) / {M^^Mp) . The Lagrangian would contain 
a term d^s^{u'^i')vu, upon insertion of the VEV of Hu- The scalars can be converted to d 
and s quarks via a gluino loop, and would result in the following estimate for processes 
such as n — )■ uK^: 

T{n -> uKy^^^^ \Yu^<i^.^Y,^^^,Y^-pHf ,^ / M' f?^)' _i_^ . (14) 

^ ' 647r ^ ^ /2 \MsJ \nJ MlMj ^ ' 

Here Ms is a typical SUSY breaking mass scale. For Fuc^c^c = YLd'^p = 10^^, Y^^-pn^ = 1.0, 
Mp = 10^ GeV, M^ = 10^6 GeV, Mg = I TeV, one obtains r(n -^ uK^) ^ 10^^ yrs. 

As for the consistency of such an intermediate scale mass for p with gauge coupling 
unification in SUSY models, we note that p(3, 2, 1/6) forms a complete SU{b) 10-plet 
along with a (3, 1, —2/3) and a (1, 1, 1) fields. If these fields also have an intermediate scale 
mass, unification of gauge coupling would work as in the MSSM. 

It should be noted that in SUSY models, the decay n — )■ e~ K^ would be suppressed, 
since the VEV of Hu picks a neutrino field in Eq. (|3]). Discovery of n — )■ e^K^ decay would 
thus hint at a deeper non-supersymmetric dynamics. 

5 Baryogenesis at the GUT epoch 

We now proceed to the computation of the baryon asymmetry of the universe induced at 
the GUT epoch. The {B — L)-violating decays of the scalars u^S, 1, —1/3) and 77(3, 1, 2/3) 
and of the vector gauge boson Vq(3,2, 1/6) will be used to illustrate the mechanism. We 
shall see that in each case, the out of equilibrium condition can be satisfied, and that there 
is enough CP violation. These decays generate an asymmetry in {B — L), which is not 
destroyed by the effective interactions induced by the electroweak sphalerons, and would 
survive to low temperatures. This is in contrast with the induced baryon asymmetry in the 
{B — L)-preserving decays of GUT scale scalars and gauge bosons in unified models such 
as SU{5), which is washed out by the spharleron interactions. 

5.1 {B — L) asymmetry in a; — >^ pH* decay 

We begin with the {B — L)-violating decay of the scalar u}{3, 1, —1/3) which is assumed to 
have a mass of order the GUT scale. [B — L asymmetry in decays of specific heavy particles 
has recently been discussed in Ref. [25].) To be concrete, we shall work in the framework 
of non-supersymmetric 50(10), although our results would hold for SUSY 5*0(10) as well, 
with some minor modifications. We identify u to be the hghtest of the various uji{3, 1, —1/3) 
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scalar fields in the 5*0(10) theory. The Yukawa couplings of Eq. ([6])- ([8]) imply that u (which 
is in general a linear combination of oJi^ul* from the I^h and 126// and 120// fields) has 
two-body decays into fermions of the type u ^ QQ, u'^e'^, u'^d'^, u^ d^, Q L. These decays 
preserve {B — L), as can be seen by assigning [B — L){u)) = —2/3. Now, u also has a 
two-body scalar decay, u — )■ pH* as shown in Fig. [6] (a), which uses the {B — L) breaking 
VEV of A'^. The scalar field p has two-body fermionic decays of the type p ^ Ld^, u'^ Q 
(the latter if kinematically allowed), which define {B — L) charge of p to be +4/3. Thus 
the decay u — )■ pH* would violate {B — L) by —2 (recall that H has zero {B — L) charge). 
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Figure 6: Tree-level diagram and one-loop corrections responsible for generating [B — L) 
asymmetry in u decay. 

Focussing on the {B — L)-violating decay u — > pH*, we define a, {B — L) asymmetry 
parameter es-i as follows. Let the branching ratio for uj — )■ pH* be r which produces a 
net {B — L) number of 4/3, and that for u* — )■ p*H be r, with net {B — L) = —4/3. The 
branching ratio for the two-fermion decays w — ;■ // is then (1— t) which has (B—L) = —2/3, 
and that for w* — )■ / / is (1 — f) which has (B — L) = 2/3. Thus in the decay oi a co + u* 
pair, a net {B — L) number, defined as es-L^ is induced, with 



es-L = {B-L)^ + {B-L)^, = -{r 
— 2(r — r) . 



f)--{(l-r)-(l-r)} 



(15) 



Note that an interplay between the {B — L)-conserving decays and the {B — L)-preserving 
decays of u is necessary for inducing this asymmetry. In addition, CP violation is required, 
otherwise r = f and thus e^-L = 0. Nonzero es-L also requires a loop diagram which has 
an absorptive part. All these conditions are realized in S'O(IO) models. The loop diagrams 
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for Lo — !■ pH* are shown in Fig. |6] (b)-(d), which involve the exchange of fermions. Since u 
can also decay to two on-shell fermions, these loop diagrams have absorptive parts. Note 
that Fig. |6] (b)-(c) are proportional to the Majorana masses for the v'^ fields, while Fig|6] 
(d) is not. These diagrams have CP violating phases, ensuring a nonzero es-L- 

We evaluate Fig. |6]in a basis where the Majorana mass matrix of the v"^ fields is diagonal 
and real. The contribution of Fig. |6] (b) to tB-L is found to be 
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(16) 



Here we have defined the trilinear scalar vertex of Fig. IH] (a) to have a coefficient \vr in 
the Lagrangian. Yq^c-p is the Yukawa coupling matrix corresponding to the coupling Q z/^ p, 
etc. Myc is the diagonal and real mass matrix of v'^ fields. Br stands for the branching ratio 
Br(a; — )■ pH*). A factor of 2 has been included here for the two SU{2)l final states in the 
decay. The function Fi{Mj) is defined as 
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with Mj denoting the mass of z/? Here G stands for the step function, signalling additional 
ways of cutting the diagram when Mj < Mp in Fig. [6] (b). 
Fig. E] (c) yields the following contribution to e^-L^ 
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where F2{Mj) is defined as 
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Fig. [6] (d) arises because in any realistic 5*0(10) model there are at least two u fields. 
The heavier u field is denoted as u'. In principle one can sum over all such u' contributions, 
but here we have kept only one such u' field. Its contribution to es-i is found to be 
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where F^{Mj) is defined as 

F3(M,) = 



Ml-M] \ (Ml 



Ml 
M2, 



(21) 



This contribution, which is non- vanishing even in the limit of vanishing Mj, requires Mj < 
M^^. Here we have defined the trilinear coupling p*u'H to have a coefficient X'vr in the 



Lagrangian. We have also assumed that M^ — M^ ^ F^, so that there is no resonant 
enhancement for the decay. (When u and u' are nearly degenerate in mass, such a resonant 
enhancement is possible. In this case, the expression for the decay rate will be smoothened 
by the width of these particles. Appropriate expressions in this case can be found in 
Ref. [26].) 

The branching ratio factor Br = Br(u; — ;■ pH*) appearing in Eqs. (TT6l) . ( TTSj) . ( I20l) can 
be estimated as follows. For this purpose let us assume that u is the field u from 10^ with 
Yukawa couplings as given in Eq. ([6]). The partial widths for the decays Ti{uj — )■ pH*) and 
^2{p — ^ //) are then given by 



Fi a; — ^ pH = 1 ^ , Fo a; -^ f f) = 



M^ . (22) 



In the expression for F2 we have assumed that u'^ is much lighter than u. In terms of 
these partial widths, the branching ratio that appears in Eqs. f lT6|) . flT6|) . fl20|) is given as 
Br = Fi/(Fi +F2). To get a feeling for numbers, let us choose a realistic set of parameters: 
M^ = 10^^ GeV, /iss = 0.6 (corresponding to the top quark Yukawa coupling at GUT 
scale) with other hij negligible, and Xv^ = (10^^, 10^^, 10^^) GeV. This would correspond 
to Br = (1.4 X 10"'^, 1.4 x 10~^, 0.58), which shows a strong dependence on Xvr. 
The total {B — L) asymmetry in u; — > pH* and its conjugate decay is given by 

^B-L - (^B^L + ("B-L + ^B-L ■ l^<^) 

This will result in the baryon to entropy ratio Yb given by 

Yb ^ "^-"^ = '-^d , (24) 

s g* 

where g^, is the total number of relativistic degrees of freedom at the epoch when these decays 
occur. In our present example g^ = 130 which includes the SM particles and the p and u 
scalar fields. The factor d in Eq. (l24l) is the dilution factor which takes into account back 
reactions that would partially wash out the induced baryon asymmetry, d is determined 
by solving the Boltzmann equations numerically, but simple analytic approximations are 
available suitable to the present setup. Defining a ratio 

T{uj ^ pH* 
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2H 
where H is the Hubble expansion rate, 



(25) 



H=1.66(7y2^, (26) 

Mpi 



the dilution factor can be written as 



d^{ ' 0.3 1^>>\\ (27) 
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These approximations work well for K < 100 or so, beyond which d would be exponentially 
suppressed. For M^ = lO^^ GeV, Xvr = (10^^ 10l^ lO^^) GeV, we find K = (1.3 x 
10~^, 1.3 X 10"^, 1.23), with the corresponding dilution factors being d = (1.0, 1.0, 0.63). 
Thus we see that there is not much dilution with this choice of parameters, although 
Br = (1.4 X 10^^, 1.4 X 10~^, 0.58) can become small for smaller values of Xvr. If we 
choose M^ = 10^^ GeV instead, and vary Xvr = (10l^ 10^^ 10^*^) GeV, then we find K = 
(0.12, 12.3, 1230) and the corresponding dilution factors to be rf = (1.0, 1.4 x 10"^, 7.5 x 
10-5), with Br = (1.3 x 10"^ 0.58, 1.0). 

Although the electroweak sphaleron interactions would not wash away the GUT scale 
induced {B — L) asymmetry, partial wash-out can occur via the {B — L)-violating in- 
teractions of the right-handed neutrinos. This is possible because the u'^ fields acquire 
{B — L)-violating Majorana masses, and their interactions with the Higgs field and the 
lepton fields can erase part of the GUT-induced asymmetry. In some cases it may be de- 
sirable to have partial wash-out. One can also prevent any wash-out by decoupling the u'^ 
fields at the same temperature as the u field (which would happen if M^c ~ M^). For the 
surviving light v^ fields, the condition Y'^ /{8n) <^ l.QQgJ {M^c/M-pi) would guarantee no 
further wash-out, where Y is the Dirac Yukawa coupling of the light v^ field. 

While we do not present the calculation of Yb in the supersymmetric version of the 
5*0(10) model, results in that case would be similar to the non-SUSY case. More diagrams 
contribute to the generation of {B — L) asymmetry with superparticle decays included. The 
number of relativistic degrees of freedom g^ would also double in this case. 

5.2 {B — L) asymmetry in 77 ^^ pH decay 

The {B — L) asymmetry induced in the decay 7/(3, 1, 2/3) — > pH is similar to the one 
induced in the decay u — )■ pH*. The tree-level {B — L)-violating decay and the one-loop 
corrections are shown in Fig. [71 rj has a fermionic decay mode t] ^ u'^u^ shown in Fig. [7] 
(a), which can be used to define its {B — L) quantum number as {B — L){i]) = —2/3. The 
decay r] — )• Hp (Fig. [11(b)) would then violate {B — L) by +2 units, since {B — L){p) = 4/3 
obtained from the decays p ^ Ld'^, u^Q. The one-loop correction to this decay is shown in 
Fig. [3 (c), which is evaluated in analogy to Fig. [6l (c) to be 



4^1l = ^Im 



Tr{F„c,c^ fJ.q^ Y^cqp M,c F^iM,.)} Xvr 



Br. (2J 



\Xvr\' 
Here the function F2{Mj) is defined as 

^.(M,)-l"(l + 4)+e(l-^)(l-^), (29) 

and Xvr identified as the coefficient of the trilinear vertex p*riH*. 
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Figure 7: Tree-level and one-loop level diagrams responsible for inducing [B — L) asym- 
metry in ?7 — 7- pH* decay. 

A noteworthy feature of baryon asymmetry generated in the decay rj — ;■ pH* is that both 
the scalars rj and p can be relatively light, since they do not induce d = 6 baryon number 
violating operators of Eq. ([T]). In supersymmetric extensions of 5*0(10), this allows for the 
possibility that the gravitino abundance problem of supergravity models can be evaded. 
Typically in supergravity models, the reheat temperature after inflation is required to be 
^reheat < 10® GcV, in Order to sufficiently dilute the gravitino abundance in the universe. 
If Mrj, Mp < 10® GeV, this requirement would be compatible with the {B — L) asymmetry 
generation. 

5.3 {B — L) asymmetry in vector gauge boson decay Vq — > V*^ H* 

The d = 7 baryon number violating operators of Eq. ([2]) can arise by integrating out the 
Vq(3, 2, 1/6) and Kjc(3, 1, —2/3) gauge bosons of 5*0(10) (see Fig. [3]), which lie outside of 
the 5*f/(5) subgroup. The decay Vq — )■ V*cH* and the conjugate decay Vq — )■ Vu^H can 
produce a primordial {B — L) asymmetry at the GUT scale, which would survive down to 
low temperatures without being washed out by the sphaleron interactions. The tree-level 
decay diagram and the one-loop correction are shown in Fig. [81 The vector gauge boson Vq 
has two-fermion decays into the following channels: Vq — )■ u'^u, Qd'^, u^u for the charge 
2/3 component, and Vq — > ud^^ u'^e, u'^d for the charge —1/3 component. These decays 
conserve (B — L), as can be seen by assigning {B — L){Vq) = —2/3. The gauge boson V*c 
has the fermionic decays V*c — )■ ed, d'^e^, Vu, u'^u'^, suggesting that {B — L)(y*c) = 4/3. The 
decay Vq — > V*cH* would then change {B — L) by +2, as in the case of the scalar decay 
uj — 7- pH* . The {B — L) asymmetry arising from Fig. [8] is found to be 



~B-L 



Im 

27r 



Tt{M,.F^{M,c)Y,clh} 



c*My 



Br, 



(30) 



with the function Fa defined as 



F,{Mj 



In 1 + 






+ e 1 




^i 






(31) 
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Here we have defined the Lagrangian coefficient of the VqVu<:H vertex to be gcMy^c- This 
is a consistent definition, since the hghter gauge boson Vuc would acquire a mass of order 
gvpt- Here c is a Clebsch factor of order unity, with its value depending on the Higgs 
representation used for rank reduction {126h or 16h)- 

^AAAAA/yVVWW 

W y. ^ U 

I * 

(a) (b) 

Figure 8: Tree graph and one-loop graph responsible for (B — L) asymmetry in the decay 
of vector gauge boson Vq. 

The factor Br that appears in Eq. f l30p is the branching ratio Br(VQ — > V*cH*). It is 
determined in terms of the two partial widths as Br = ri/(ri + r2), where 

T2iVQ^ff) = ^Mv^. (32) 

If we choose Mvq = 10^^ GeV, K = Ti/H ~ 0.016 (with g = 0.55 and c = 1), so that the 
dilution factor is d ^ 1.0. There is a modest suppression in e^_^ arising from the branching 
ratio, since Br ~ 0.06 with this choice. 

5.4 Baryon asymmetry in a class of minimal 50(10) models 

We now show how the GUT scale induced asymmetry in cj — )■ pH* decay can consistently 
explain the observed value of Yb = (8.75 ± 0.23) x 10~^^, in a class of minimal 5*0(10) 
models. In this class of models, a single 10^ and a single 126^ couple to fermions, as in 



Eqs. (I6])-(I7]). In SUSY models this is automatic with a single IOh and 12Qh employed. 
In non-SUSY models the 10^ can also couple (10^:^ must be complexified to generate real- 
istic fermion masses), however if a Peccei-Quinn symmetry is assumed, the 10|^ coupling 
would be absent. This class of models is highly constrained due to the small number of 
parameters that describe the fermion masses and mixings, and leads to predictions for the 
neutrino oscillation parameters p]. In addition to generating large mixing angles for solar 
neutrino oscillations and for atmospheric neutrino oscillations, these models predict a rela- 
tively large value of ^13, viz., sin^ 26'i3 ^ (0.085 — 0.095), both in the non-supersymmetric 
and the supersymmetric versions, consistent with recent results from Daya Bay and other 
experiments [S]. 
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To illustrate how realistic choice of parameters generate acceptable Yb, we choose the 
u field to be almost entirely in the IOh- We also choose X'vr that appears in Fig. [6] (d) to 
be small, so that the leading contribution to es-L is from Fig. |6] (c), as given in Eq. f lTSj) . 
In this limit, we find 



2V3\h3f: 



|2 




e^-L^^M^U + ln 1 + T7f sin0. (33) 



Here we have kept only the third family Yukawa couplings, which is the leading contribution, 
and we have defined (p = a.rg{hl^f^X + |}. Choosing h^^ ^ 0.6 (the top quark Yukawa 
coupling at the GUT scale), and A = 0.25, vr = 10^^ GeV, /g = lO^^ (so that f-^VR = 10^^ 
GeV, consistent with the light u^ mass arising via the seesaw mechanism), = 0.12, we 
find eB-L = 1-6 x 1.9 x 10"^. If M^ = 10^^ GeV, then Br = 0.96, K = 197 so that the 
dilution factor is d = 5.6 x 10~^. This results in a net Yb = 8.2 x 10"^^, consistent with 
observations. 

While natural choices of parameters can generate acceptable Yb, due to the high sensi- 
tivity of dilution factor on the masses of the heavy particles, precise predictions are difficult 
to make. For Vq — )■ V*cH* decay we find the process to be typically out of equilibrium so 
that (i ^ 1 for My^ ~ 10^^ GeV. Natural values of the asymmetry parameter in this case is 
eB-L ~ 10^^. Some dilution effects from the z/'^ interactions would be welcome in this case. 

It should be mentioned that the d = 7 operators of Eq. ([5]) also arise naturally in 
S'O(IO) models, as already noted. The r]*pH and the ^pH* vertices arising from the 126^ 
couplings can be used for GUT scale {B — L)-genesis without generating d = 7 nucleon 
decay operators. The decays rj — > pH has already been analyzed, but if these particles arise 
from 126h they do not lead to {B — L)-violating nucleon decay. 

6 Conclusion 

In conclusion, we have pointed out that the complete set oi d = 7 baryon number violating 
operators that lead to the selection rule A{B — L) = ±2 in nucleon decay can emerge as 
effective low energy operators in S'O(IO) unified theories with either a 126// or a 16// Higgs 
field used for breaking the B — L gauge symmetry. The strength of these operators is unob- 
servable in single-step models where 5*0(10) breaks directly down to the standard model. 
In non-supersymmetric 5*0(10) models, an intermediate symmetry is required in order for 
the gauge couplings to unify correctly. We have shown that in several instances with such 
an intermediate scale, the d = 7 baryon number violating operators can lead to observable 
nucleon decay rates. The decay modes are distinct from the conventional GUT-motivated 
modes, and include n — )■ e~K~^, e~7r+, etc. We have also identified supersymmetric scenar- 
ios where such modes may be within reach of experiments, consistent with gauge coupling 
unification. 
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A second major result of this paper is a new way of generating {B — L) asymmetry 
in the early universe by the decay of GUT mass particles. It is these particles which also 
induce the d = 7 nucleon decay operators. Such an asymmetry is sphaleron-proof, in 
that it does not get erased by the effective interactions of the electroweak sphalerons. We 
present several examples where consistent asymmetry can be generated with the GUT scale 
decays of particles obeying the A{B — L) = ±2 selection rule. Further, we show that in 
minimal 5*0(10) models which explain the large neutrino mixing angles and predict relative 
large value for 6*13, consistent with recent experimental results, that the induced baryon 
asymmetry via the proposed GUT-scale mechanism is compatible with observations. There 
is thus a strong connection between neutrino oscillation parameters and baryon asymmetry 
in this class of models. 
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